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ABSTRACT 

We make the hypothesis that molecular cloud fragments are triaxial bodies with a large scale magnetic 
field oriented along the short axis. While consistent with theoretical expectations, this idea is supported 
by magnetic field strength data, which show strong evidence for fiattening along the direction of the 
mean magnetic field. It is also consistent with early submillimeter polarization data, which show that 
the projected direction of the magnetic field is often slightly misaligned with the projected minor axis 
of a molecular cloud core, i.e., the offset angle \E' is nonzero. We calculate distributions of for various 
triaxial bodies, when viewed from a random set of viewing angles. The highest viewing probability 
always corresponds to ^' = 0°, but there is a finite probability of viewing all nonzero 5*, including even 
^ = 90°; the average offset typically falls in the range 10° — 30° for triaxial bodies most likely to satisfy 
observational and theoretical constraints. 

Subject headings: ISM: clouds - ISM: magnetic fields - MHD - polarization 



1. INTRODUCTION 

Observations establish that magnetic fields play an im- 
portant role in shaping the structure and dynamics of 
molecular clouds. Polarization maps reveal that there is 
a magnetic field threading molecular clouds that is or- 
dered over large scales (e.g., Vrba, Strom, & Strom 1976; 
Goodman et al. 1990; Schleuning 1998) and Zeeman mea- 
surements establish that measured field strengths are in- 
variably strong enough to yield a mass-to-flux ratio very 
close to the critical value (Crutcher 1999; hereafter C99), 
i.e., the magnetic energy approximately equals the gravita- 
tional binding energy. A problem with absorption maps of 
the polarization of background starlight is that the dens- 
est regions (the cloud cores) are not sampled (Goodman 
et al. 1995). However, new maps of polarized submillime- 
ter emission (Matthews & Wilson 2000; Ward-Thompson 
et al. 2000) do sample these regions and provide a cru- 
cial test of theoretical scenarios involving magnetic fields; 
in particular, the relative alignment of the projected mag- 
netic field with the projected minor axis of the core density 
distribution is an important diagnostic. 

When making the above comparison, two important 
factors must be kept in mind: 1) the measured relative 
alignment is only between two projected quantities in the 
plane of the sky, and care must be taken to draw infer- 
ences about the relative orientation in the actual three 
dimensional object, and 2) theoretical models of magneti- 
cally dominated clouds (e.g., Mouschovias 1976; Lizano & 
Shu 1989; Fiedler & Mouschovias 1993), while predicting 
the important effect of flattening along the mean mag- 
netic field direction, are usually restricted (for numerical 
convenience) to the physically unnecessary assumption of 
axisymmetry. In fact, it has been known for a long time 
(Mestel 1965; Lin, Mestel, & Shu 1965) that gravity am- 
plifies anisotropic structure, so that an initially unstable 
fragment will first collapse along one dimension, forming 
a sheet, which will subsequently break into elongated fila- 
mentary structures. The presence of a large scale magnetic 
field can enhance this effect, as the first stage of collapse 



is preferentially along the mean field direction. The sub- 
sequent contraction lateral to the field will occur in qual- 
itatively the same manner as described by Mestel (1965), 
particularly if the cloud is magnetically supercritical, but 
even if it is subcritical, since the evolution is still driven by 
gravity (Mouschovias 1978; Langer 1978), but on the am- 
bipolar diffusion (rather than free-fall) time scale. Thus, 
the most general formation mechanism for a dense core 
in the magnetic field scenario yields an object flattened 
along the mean fleld direction and also having an elon- 
gated structure in the lateral direction; such an object is 
triaxial. 

The properties of the projected shapes of triaxial bod- 
ies have been investigated in the context of galaxy studies. 
Stark (1977) has proven the important result that the iso- 
contours of a triaxial body seen in projection are elliptical. 
Binney (1985) has extended this analysis, and his results 
can be applied to find the distribution of angular offsets 
5* between the projected field direction and the projected 
minor axis of a triaxial body, given that the body is pri- 
marily fiattened along the direction of the mean magnetic 
field. 

In § 2, we review observational evidence to establish that 
molecular cloud fragments are ubiquitously flattened along 
the mean magnetic field direction, and in § 3 calculate the 
distribution of ^ for various triaxial bodies. In § 4 we dis- 
cuss how the results relate to current observations and to 
other theoretical scenarios presented in the literature. 

2. CLOUD FLATTENING ALONG MAGNETIC FIELD LINES 

The magnetic field strength data of C99 also reveals 
that the correlation of the line-of-sight field strength i3ios 
with the density p is in apparent agreement with models of 
preferential flattening along the magnetic fleld {B cx p^^^] 
Mouschovias 1976) rather than isotropic contraction {B cx 
p^/^). Figure la (equivalent to Fig. 1 of C99) plots (in 
logarithmic form) the conflrmed magnetic fleld strength 
detections of C99 against the measured number densities 
n. The solid line is the least squares best fit, with slope 
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0.47 ±0.08, and chi-squared error statistic = 2.36. The 
correlation coefficient is 0.84. 

We seek an even better correlation for the C99 data as 
follows. Consider a cloud that is flattened along the mean 
magnetic field direction. In this direction, force balance 
between an isotropic pressure and self-gravity yields 



(1) 



where po is the density at the midplane, Uv is the total 
(thermal and non-thermal) one-dimensional velocity dis- 
persion, S is the column density of the cloud, and G is the 
gravitational constant. Here we assume that the effect of 
an external confining pressure is negligible. Furthermore, 
the relation 



|^p(2.Gi/2)-i 



(2) 



defines the dimensionless ratio /i of the mass-to-flux ra- 
tio to the critical value {2'kG^/'^)~'^ for a disk (Nakano & 
Nakamura 1978). Combining equations (Q) and (Q) yields 



B = (87rci) 



1/2 ^vP 



1/2 



(3) 



(see Mouschovias 1991 for a similar relation), where ci (^ 
1) is an undetermined proportionality constant between 
Po and the mean density p. Clearly, the correlation in 
Figure la is imperfect because variations in a^, and fj, are 
unaccounted for, and the observed quantity B\os = B cos 6 
(where 9 is the angle between the field direction and the 
line of sight) does not correspond exactly to B. We im- 
prove the correlation by accounting for variations in the 
observable quantity Uy. Figure lb plots -Bios versus n^/^ 
in logarithmic form; equation (^) predicts a slope of 1 to 
the extent that p, is constant from cloud to cloud. We 
use values of A^; = cr„(81n2)i/2 listed in Table 1 of C99. 
The correlation coefficient in this plot is 0.95. The solid 
line is the best least squares fit and has slope 1.00 ± 0.09, 
and — 0.79. The dashed line is the theoretical relation 
(|) for yjoil p- = 1. The two fines coincide if yjclj \i — 0.8. 
This exceptional correlation leads us to the conclusion that 
clouds (and cloud fragments) are flattened (at least mod- 
erately) along the magnetic field direction, even in the 
presence of turbulent motions, due to the dynamically im- 
portant roles of self-gravity and the mean magnetic field. 
We note that the correlation found in Figure lb is equiva- 
lent to that presented by Myers & Goodman (1988) in their 
Figure 1 (see also Mouschovias & Psaltis 1995) based on 
earlier magnetic field data. We also note that our fit to the 
magnetic field data is quite general; it does not require ax- 
isymmetry, or that the Alfven speed va = B / be 
constant from one cloud to another. However, equation (|3|) 
implies that a^jvA = (2ci)~"'^/^/z, so that the Alfven mach 
number may be approximately constant, if non-thermal 
motions dominate thermal motions and /i does not vary 
much from one cloud to another. 
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Fig. 1 a) Plot of log -Bios versus log n from the C99 data 
set. The solid line is the best fit (see text). The correla- 
tion coefficient is 0.84. b) Plot of log Bios versus Xoga^v}!"^ 
for the same data. The correlation coefficient is 0.95. The 
solid line is the best fit, and the dashed line is a theoretical 
relation (eq. H) for ^/cl/ p. = 1). 

3. PROJECTION OF TRIAXIAL BODIES 

Let z be the direction of the mean magnetic field. Since 
the mean field contributes cloud support against gravity 
only in the lateral [x and y) directions, but the cloud need 
not be axisymmetric in the x — y plane, the most general 
cloud configuration is a triaxial body with the shortest 
dimension along the z— axis. The simplest mathematical 
description of such a body is 



2 , y 2 



(4) 



where a is a constant, and !>(,>£,■ Here we follow the 
geometric analysis of Binney (1985) to calculate projected 
quantities when the object is viewed from an observing 
angle (^, (j)) (where the angles are defined on an imaginary 
viewing sphere and have their usual meaning in a spherical 
coordinate system). Defining the quantities 



A = 



cos 



COS 



sm 
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B = cos 6* sin 2(/) ( 1 - — J — , 



1 \ 1 



and 



C . 



1 



(5) 
(6) 

(7) 
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it can be shown that the projection in the sky of the tri- 
axial body is an elhptical body of apparent axial ratio 



A + C - D 
A + C + D 



1/2 



(8) 



where D = ^/{A — Cy + B"^. Furthermore, the quantity 
tp = \ arctan ( —r——^] (9) 



can be used to find the absohitc vahic of the angular off- 
set between the projected z axis and the apparent minor 
axis of the projected ellipse. Letting p= (A — C) cos 2^ + 
Bsin2^, the required positive definite angle is 



1^1 
7r/2- 



for p < 
for p > 0. 



(10) 



Since 'i> corresponds to the offset between the projected 
magnetic field B± and the projected minor axis of a core, 
we calculate the distribution of observed offsets using a 
Monte Carlo simulation. A large number N = 10^ of ran- 
domly oriented viewing angles (0, (p) are used to generate 
a distribution of 4* for a given set of intrinsic axial ratios 
C and ^ of a triaxial body. Figure 2 shows, for each of 
three cases ^ = 0.1,0.3, and 0.5, the probability distribu- 
tion function / versus ^ for three different values of (, as 
well as the dependence of the mean value (4') of the distri- 
bution and the mean axial ratio (q) on ^. The lower limit 
of C = ^ corresponds to a prolate object and the upper 
limit ( = 0.9 is very nearly an oblate object. 

Our calculated distributions show the robust result that 
there is always a peak in / at ^' = 0, so that it is the most 
probable value, but that there is a long tail towards all 
nonzero angles, including ^ = 90°. Figures 2b, 2d, and 
2f show how (4') varies for clouds of different triaxial ra- 
tios; it becomes smaller as a cloud goes toward the oblate 
limit. For a perfectly oblate cloud (C = 1), there would 
be no distribution, as all viewing angles yield = 0. The 
distributions for prolate clouds (( = have the largest 
possible values for (^'). The values of (q) illustrate which 
triaxial models may be most likely. The crosses highlight 
the values of ( for which (q) is in the range 0.5 — 0.65 that 
is consistent with observations of dense cores (e.g., Myers 
et al. 1991; Jijina, Myers, & Adams 1999). 

While data sets that measure q and for a large num- 
ber of cores will ultimately settle the values (and distribu- 
tions) of C and ^, we do not currently favor highly flattened 
cases ^ w 0.1 due to 1) the observational result that highly 
elongated cores arc rarely observed, and 2) the theoreti- 
cal result that flattening along the magnetic fleld yields 
cores of axial ratio « 0.25 — 0.33 (see Ciolek & Basu 2000 
and references therein). Hence, the ^ = 0.3 curves repre- 
sent likely scenarios, although the ^ = 0.5 curves cannot 



be ruled out observationally. The observational constraint 
on (q) also allows the; widest plausible range of values 
when ^ = 0.3, rather than imply near-oblate (for ^ = 0.1) 
or near-prolate (for ^ = 0.5) configurations. For ^ in the 
range 0.3 — 0.5, clouds with (q) in the range 0.5 — 0.65 have 
(*) in the range ~ 10° - 30°. 

To illustrate the point that no single observation can es- 
tablish the relative orientation of the fleld and the minor 
axis of the intrinsic three-dimensional object, we present 
in Figure 3 a view of a triaxial body with (^, C) = (0.3, 0.6) 
from three positions on the viewing sphere. This illustrates 
that all possible orientations of B± with the projected mi- 
nor axis are possible, although the first panel shows the 
single most likely orientation and the middle panel is clos- 
est to the average orientation (^) = 20° for this case. 

4. DISCUSSION 

Our predicted probability distributions, with their cor- 
relation toward 4" = 0, are in good agreement with early 
submillimeter polarimetry of dense regions. Matthews & 
Wilson (2000) measure three distinct regions in OMC-3, 
of which two imply near perfect alignment of B± with 
the projected minor axis, and one has ^ ^ 30°. Ward- 
Thompson et al. (2000) infer B± in three dense cores, 
LI 544, LI 83, and L43, and also find it to be correlated with 
the projected minor axes. They measure ^ = 29°, 34°, and 
44° for the three cores, respectively, although the last mea- 
surement is complicated by a nearby outflow and a weak 
polarization signal toward the core center. For L1544, the 
nonzero \1/ implies that the ax;isymmetric assumption in 
the magnetic model for that core developed by Ciolek & 
Basu (2000) needs to be relaxed, at least moderately. How- 
ever, since that model does already account for flattening 
along the field lines, the magnetic field strength predic- 
tions would likely not change substantially. 

Future large data sets of the offset angle ^P, as submil- 
limeter polarization measurements become common, will 
properly test the distributions of 5" presented in this paper. 
Such data will provide a clear observational test between 
the scenario presented here (which is already supported 
by the Zeeman data) , in which gravity and the mean mag- 
netic field play an important role in shaping core struc- 
ture (even in the presence of turbulence) and those sce- 
narios in which turbulent compression is dominant, lead- 
ing to the claim (Ballcsteros-Parcdes, Vazquez-Semadeni, 
& Scalo 1999) that there is no correlation between the 
magnetic field direction and the cloud elongation. 



I thank the anonymous referee for comments which 
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by a grant from the Natural Sciences and Engineering Re- 
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^Such clouds axe distinct from existing models of prolate clouds (Tomisalsa 1991; Fiege & Pudritz 2000) in which the mean magnetic field 
lies along the major axis; here the mean field lies along the minor axis. The cloud can be prolate (C = by chance, but there is no intrinsic 
symmetry that favors such a configuration, and the pressure exerted by the mean field makes it likely that C > 4 for most clouds. 
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Fig. 2 Probability distribution function / versus angle ^l* for various values of (^,C)j well as the dependence of mean 
values (^t') and (q) on ( for individual values of ^. a) / versus for ^ = 0.1 and ( — 0.1,0.5,0.9. b) (4') (solid line) and 
(g) (dashed line) versus ^ for ^ = 0.1. c) Same as a) but ^ = 0.3 and ( = 0.3, 0.6, 0.9. d) Same as b) but ^ = 0.3. e) Same 
as a) but ^ = 0.5 and C = 0.5,0.7,0.9. f) Same as b) but ^ = 0.5. In b), d) and f), the crosses highhght the region in 
which (q) = 0.5 — 0.65, consistent with observations. 
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Fig. 3 Simulated contours (solid linos) and mean polarization direction (dashed lines) for a triaxial body with axial ratios 
= (0-3,0.6) seen from three sets of viewing angles {0,(j)). The The x' and y' axes are in the plane of the sky and 
are chosen to lie along the projected major and minor axes, and a' is the apparent semi-major axis. The dashed lines lie 
along the projected direction of the intrinsic z axis, a) {9, (p) = {n/2,0), yielding an apparent axis ratio q = 0.3 and offset 
angle * = 0°. b) {9, <j)) = {tt/A, 7r/4), yielding q = 0.57 and ^ = 28°. c) {9, 0) = (tt/G, 0), yielding q = 0.68 and * = 90°. 



